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Abstract 



o 

(N 

• We prove that in any static spacetime the maximization of matter entropy subject to the continuity 

D I equation for the stress-energy tensor and the initial- value constraint leads to the full set of Einstein 

equations and to Tolman's relation between local temperature and lapse function. This result is shown 
00 , to imply that the thermodynamics of gravitating systems in equilibrium is fully specified by variables 

defined on the system's boundary, in particular, the boundary's geometry and extrinsic curvature. 
Hence, self-gravitating systems in equilibrium are holographic already at the level of classical relativity. 

O^' The association of entropy to black holes by Bekenstein and Hawking strongly suggests of a fundamen- 

^ ' tal relationship between gravity and thermodynamics. It has led to the conjecture that the fundamental 
theory of gravity should satisfy the holographic principle [l] , namely, the statement that the full informa- 
tion about a gravitational system is contained in the degrees of freedom of the system's boundary. 

In this letter, we show that self-gravitating systems in general relativity are holographic already at 
. the level of the classical theory with no reference to black holes or the Bekenstein-Hawking entropy. We 
' demonstrate that the thermodynamical properties of a self-gravitating system in equilibrium are defined 
solely in terms of variables defined on the system's boundary, in particular, the boundary's metric and 
_ ■ extrinsic curvature. 

, To this end, we first prove that in a general static spacetime the full set of Einstein's equations follow 

' from the maximization of matter entropy subject to (i) the continuity equation for the stress-energy 
tensor, (ii) the initial value (Hamiltonian) constraint and (iii) suitable boundary conditions. This proof 
generalizes previous results that were derived for spherically symmetric static spacetimes — see, Refs. 
^ I [21 E] for specific model systems, Ref. [4] for the most general spherically symmetric case, and Ref. [5] for 
H ' elaborations. 

The proof proceeds as follows. Consider a static globally hyperbolic spacetime M = R x T, with 
four-metric 

ds'^ = -L^{x)dt^ + hij{x)dx^dx^ , (1) 

expressed in terms of the spatial coordinates x* and the time coordinate t. The time-like vector field 
= [d/dt)^ is a Killing vector of the metric Eq. ([1]). L is the lapse function, and hij is a t-independent 
Riemannian three-metric on the surfaces of constant t. The time-like unit normal on St is = Ld^t 
and the extrinsic curvature tensor on vanishes. 

Thermodynamic definitions. Let V C S be a compact spatial region, with boundary B = dV . V 
contains an isotropic fiuid in thermal and dynamical equilibrium, described by the stress-energy tensor 

T^u = p{x)n^ny + P{x){g^y+n^n^), (2) 

where p{x) and P{x) are the energy density and the pressure, respectively. 



* anastop@physics.upatras.gr 
t ksavvidouQphysics . upatras.gr 



1 



The stress-energy tensor satisfies the continuity equation V^T^'^ = which implies that 

V,;P V,;L 



p + P L 



(3) 



Eq. ([2]) is the standard form of the stress-energy tensor for ideal fluids. Since dissipative processes 
are absent in equilibrium configurations, Eq. ([2]) also applies to non-ideal fluids. 

We assume that the fluid consists of r different particle species. The associated particle-number 
densities na{x), a = 1, . . . , r, together with the energy density p{x) define the thermodynamic state space. 
All local thermodynamic properties of the fluid are encoded in the entropy-density functional s{p,na). 
The first law of thermodynamics takes the form 

Tds = dp - ^ Padua, (4) 

a 

where pa = is the chemical potential associated to particle species a and T = ^f^^ is the local 

temperature. The pressure P is defined through the Euler equation 

p + P-Ts-^/i^na = 0. (5) 

a 

Combining Eqs. ([5]) and (jH), we derive the Gibbs-Duhem relation, dP = sdT + ^^Uadpa- 

Entropy maximization. Next, we maximize the total entropy S = jy d^x\/Tis{p, Ua) for fixed values 
of the particle numbers Na = Jy d^x\/Tina{p,na)- Entropy maximization is subject to the continuity 
equation Q for the fluid, and to the Hamiltonian constraint of general relativity. For static spacetimes, 
the Hamiltonian constraint reads 

C{x) := 16ttp{x) - R{x) = 0, (6) 

where R is the Ricci scalar associated to the three- metric hij. 

The Hamiltonian constraint expresses the energy density /j as a function of the metric hij and the 
continuity equation expresses the lapse L in terms of the energy density p. It follows that the entropy 
S and the particle numbers Na are functionals of the three-metric hij and the particle-number densities 
Ua- Entropy maximization requires the variation of + Yla^<i^a, with respect to Ua and hij for some 
Lag range multipliers b^. 

Variation of 5 + J2a ^aNa with respect to yields. 



^s = Y.f ^^^^ + ^'^- = nf ^"^^^ (- Y + ^'^'^ " °' 



(7) 



leading to 5^ = Hence, for equilibrium conflgurations the thermodynamic variables ^ are constant 
in V. 

We introduce the function uj (a Massieu function [7J) as the Legendre transform of the entropy density 
s with respect to Ua 

X - ds \ p + P . 

uj{p,ba) := s - 2_^- na = S + } ^ baUg = ^ . (8) 

a a 

Expressed in terms of oj, Eq. @ takes the form 

duj = Y -"^^adba- (9) 

a 

It follows that = d(jj/dp and Ua = —duj/dba- The Gibbs-Duhem relation becomes dP = iodT + 
TY^gUadba- 
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Tolman's relation. Since the variables ba are constant for entropy-maximizing configurations, dP/dT = 
uj = {P + p)/T in V. Combining with Eq. ([3]), we obtain 

— = -— • <i°> 

which leads to Tolman's relation between local temperature and lapse function [6] 

LT = n, (11) 

where T* is a constant. In an asymptotically flat spacetime, L = 1 at spacelike infinity. Then, T^, is 
identified with the temperature as seen by an observer at infinity. Eq. (jlip is the generalization of the 
zero-th law of thermodynamics: the temperature of a self-gravitating body in equilibrium is everywhere 
constant when transformed to the frame of an observer located at infinity. We note that the only as- 
sumptions required for the derivation of Tolman's relation are (i) that the entropy density s does not 
depend on the metric hij and (ii) that the Hamiltonian constraint does not involve the particle densities 
Ua- Tolman's relation is essentially kinematical, in the sense that its derivation does not depend on the 
explicit form of Einstein's equations. It follows from Galileo's principle that gravity is insensitive to the 
physical composition of material bodies. 

Derivation of Einstein's equations. By employing the definition Eq. ([8]), entropy-maximization is 
equivalent to the maximization of the functional fi, defined as 

n[hij,ba] := j d^xVhu}{p,ba), (12) 

The variation of 0, with respect to metrics hij that satisfy the Hamiltonian constraint Eq. Q yields 



6n 



J d'xVh (^^h^^h, + ^^p) = I d'xVh (^^h^^6h,, + ^) . (13) 



Using the equation 5R = —W'^Shij + V* {V^dhij — h^^^iShki) for the variation of the Ricci scalar R, 
together with Eq. PU]) . we find 



(14) 



Sij Rij _ hijR _ 1(V*V^L - h'^VkV'^L) - Snh'^P. (15) 



where 



In order to compute the total divergence term in Eq. ()14p . we perform a 2 + 1 decomposition of the 
metric hij on the bounding surface B. Let B be described locally by the condition f{x) = 0. The unit 
normal to B is rrii = aS/if, where a = \ j \Jh}^WiJWj] . The induced metric on B is Oij = hij — miruj 
and the extrinsic curvature tensor of B is Kij := a^ajV^mi. 

Then, Eq. ([H]) becomes 



(50 = - y d^xVhS'^hij + ^ d^y^/^ [m}ViLa^^5aki - 2La'^Kij + Ln^^aij^ , (16) 

where the coordinates on B are denoted as y. 

For a variation with fixed values of the two-metric dij and of the extrinsic curvature Kij on the 
boundary, 5^1 = — J d^xVhS^^ 5hij . Thus, the maximum entropy principle, dfl = 0, leads to the condition 
S"^^ = 0, which coincides the spatial components of Einstein's equations for the static spacetime metric, 
Eq. ^. 
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The equation S'*-' = corresponds to Hamilton's equations for the momentum vr*-', conjugate to hij 
[8]. In the Hamiltonian formulation of general relativity, the equations of motion incorporate the sym- 
plectic structure (Poisson bracket) and the constraints of the gravitational state space. The Hamiltonian 
constraint was employed in the derivation of Eq. ()16p . Hence, we conclude that the thermodynamic de- 
scription of the system contains implicitly information about the symplectic structure of general relativity 

m- 

The thermodynamic state space of self- gravitating systems. Next, we employ Eq. ([16]) in order to 
formulate the thermodynamic properties of self-gravitating systems. To this end, we concentrate on 
solutions to Einstein's equations, S^^ = 0, and examine the changes to under external variations of the 
boundary conditions. 

Eq. ([16]) applies to any spatial region V and boundary B. When studying an isolated self-gravitating 
system, B is to be identified with the system's physical boundary, for example, a bounding box (as in 
ordinary thermodynamics) or a stellar surface. B separates between a region in which the metric hij is 
given by a solution of Einstein's equation with matter (the interior) and a region where hij is a solution 
of the vacuum Einstein equations (the exterior). Assuming that the metric hij is regular in the interior 
region (i.e., everywhere locally Minkowskian, no inner boundaries), then Eq. (jl6p leads to 

5n = r^^-^ j (fy^/^ [m'ViLa^^Saki - 2La'^Kij + Ln'^daij^ . (17) 

We observe that Q varies only with the two-metric aij and the extrinsic curvature Kij of the boundary. No 
other components of the three metric and its derivatives have thermodynamic significance. O also depends 
on the variables b^, these are constant in V, and thus their value is also determined at the boundary. It 
follows that the thermodynamic state space for a self-gravitating system in equilibrium consists only of 
variables that are accessible at the system's boundary, and that the function Q[aij, Kij ,ba], defined solely 
on boundary data, fully specifies the thermodynamic properties of the system. 

Moreover, by Eq. (fTTll the thermodynamically conjugates 6il/6aij{s) and 6VL/5nij{s) are also boundary 
variables. The particle numbers Na = are bulk variables, however, they are conserved in a closed 

system, and their value can be ascertained from information about the process through which the system 
was formed. An external observer with no knowledge about the system, other than the conserved numbers 
Na, can fully reconstruct the equations of state (the functional relation between state-space variables and 
their conjugates) solely by local measurements at the boundary. Thus, the thermodynamics of self- 
gravitating systems exhibit a striking manifestation of the holographic principle; all hydrodynamic bulk 
variables (density, pressure) are encoded into geometric properties of the boundary. 

We must emphasize that the holographic description, derived here, refers only to thermodynamic 
properties of a self-gravitating system; not all bulk variables can be identified by boundary variables. 
This is the case in particular for the three metric hij on V. The thermodynamic variables aij and Kij 
correspond to three components and three first derivatives of hij. Viewed as boundary data, they are 
not sufficient for finding a unique solution to Einstein's equations in the interior. Hence, many different 
solutions to Einstein's equations correspond to the same values of the thermodynamic variables on the 
boundary. 

This point has been made in Ref. |10j for a particular class of spherically symmetric solutions cor- 
responding to self-gravitating radiation in a box. It was shown that, nonetheless, the maximum entropy 
principle leads to a unique solution to Einstein's equations solely from the knowledge of the boundary 
variables. In this system, most solutions of Einstein's equation with fixed boundary condition contain 
a conical singularity at the origin; regular solutions form a set of measure zero in the space of all solu- 
tions. Including a contribution from the singularity to the total entropy of the systems allows for the 
implementation of the maximum entropy principle. Regular solutions turn out to maximize entropy for 
any fixed values of the thermodynamical variables on the boundary. We believe that the methodology 
of Ref. [To] can be generalized to a larger class of spacetimes, possibly to any self-gravitating system in 
equilibrium. This requires the extension of the maximum entropy principle also to non-regular spacetimes, 
taking into account the contribution of internal boundaries and singularities. This issue is presently under 
investigation. 
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Weak gravity limit. It is instructive to see how the correspondence between bulk-hydrodynamics and 
boundary geometry works at the hmit of weak gravitational interaction. We assume that the boundary 
corresponds to a bounding box. For low densities and small system size, we expect to recover ordinary 
thermodynamics expressed in terms of extensive variables. Hence, the matter density p is approximately 
constant on the interior of the box. By Eq. the Ricci scalar R will also be constant. Hence, the 
geometry inside the box will correspond to a region of a homogeneous three-sphere of radius y^6/R. For 
simplicity, we assume a spherical boundary, so it is convenient to use coordinates {r,6,(p), in which the 
three-sphere metric reads 

ds"^ = ^ + r^{de^ + sin^ edcp'^). (18) 



At the weak gravity limit, R is close to zero, so Eq. (jlSp approaches the Minkowskian metric. For a 
boundary B defined as r = rg, aij is the standard two- metric for a sphere of radius tq. The extrinsic 

curvature is Kij ~ k^^^ ^1 — ^^3"'° )^ 1 where k-^^ is the extrinsic curvature of the surface r = tq when 
embedded in Minkowski spacetime. Thus, the energy density of ordinary, extensive thermodynamics 
is proportional to the deviation of the extrinsic curvature from its Minkowskian value: p = —^^{k — 

where k = a^^ Kij. The function ~ (jj{p, &a)(^'^o) expressed solely in terms of the geometric 
properties ro and (k — k^^^)/k^^'^ of the boundary. 

First law of thermodynamics. The projection of the equations of motion along m*, Sijm^m^ = 0, 
implies that 



m'ViL = - 

K 



SttP + yCTi/ViVjL + liKijK'^ -K^ + ^R) 



(19) 



where ^R is the curvature scalar on B and ^Vj the covariant derivative associated to oij. A natural 
thermodynamic assumption is that the local temperature is constant on i?, so that there is no boundary 
heat flow. This imphes that ^V^L = 0. Substituting Eq. ([HD into Eq. pT]) . we obtain 

If IP 

Sn = id^s—5{^)-y^NJba 
1^ J K ^ 

+ j dhL^ {n^^5a., - 2a^^5n., + ^12!l^—^±^a''^6a,i^ , (20) 



where we included variations with respect to ba- Eq. (j20p is the first law of thermodynamics for a 
self-gravitating system. The surface integral involving P corresponds to mechanical work generated by 
variations of the boundary area. The heat term T^6S and the chemical potential terms Yla l^o-^-^a are 
readily obtained from T*(5J7 + ^^a^a^^a)- However, the identification of internal energy is not obvious. 
Let us denote the term the second line of Eq. (|20p by Q/T^; is an one- form on the thermodynamic 
state space. The exterior derivative of does not vanish, so is not a exact form. Hence, it cannot be 
identified with a variation of the internal energy; it also involves work terms. 

Eq. (]20p simplifies significantly for spherically symmetric systems. The exterior metric is Schwarzschild, 
by virtue of Birkhoff 's theorem, and it is characterized by a value M of the Arnowitt-Deser-Misner (ADM) 
mass; (Jij is the two- metric for a sphere of radius rg; Kij = V^l -2M/ro/roaij, and L = v^l - 2M/ro . 
Then, = IGvrJM, and the first law of thermodynamics takes the transparent form 



5M = TJS - P(ro)(47rrg(5ro) + L{ro)Y, l^M^Na. (21) 

a 

In spherically symmetric systems, internal energy is identified with the ADM mass, and the one-form 
contains no work terms. This suggests that, in the general case, the work terms in correspond to 
deformation of boundary inhomogeneities. 
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The Komar mass. For static solutions to Einstein's equations, a natural candidate for the internal 
energy is the Komar mass M, defined in terms of the surface integral 



M = -3- / (fy^/am'^VkL. (22) 

for any closed surface B' enclosing B, since the energy density and pressure vanishes in the exterior of B. 

We will show that the Komar mass is closely related to the thermodynamic function Q. For a grav- 
itating system bound by a box, the pressure is discontinuous at B. The metric components and the 
derivatives corresponding to the thermodynamical variables aij and Kij are continuous across the bound- 
ary, but the derivative m^VkL exhibits a jump, which we will denote as [m^'VkL] . This jump corresponds 
to a non-zero value of a surface stress-energy tensor on the boundary [TI]. From Eq. (\W\i . we obtain 
[m'^Vfc-L] = — Svr-^. Taking the limit of B' ^ B from the outside, Eq. ([22]) becomes 

M = -2(f d?y^/^— + % f d^xVh ^^^^ , (23) 
Jb 1^ Jv T 

In deriving Eq. (p3]) . we employed the equation 47r(p + 3P)L = VfeV'^L, that follows from Sijh^^ = 0, 
and Tolman's law, Eq. (fTT|) . Using Eq. (fT2]l . we find 

(24) 



Jv T n 



For a fluid with a linear equation of state P = 7/>, where 7 is a constant, Jy d^x^/hp/T = 1^/(1 + 7). 
Then, Eq. (p^ becomes 

^ 1 + 7 M + 2/^(i2j/V^LP/^ 

1 + 37 ■ ^ ^ 

Eq. (j25p demonstrates explicitly how O is expressed solely in terms of variables that are defined on the 
boundary B. For spherical symmetric systems, Eq. ()25p reproduces the results of Ref. [2j. 

For a general equation of state, a closed expression for $7 in terms of the Komar mass is not possible, 
but we can identify an upper and a lower bound. For a fluid that satisfles the weak and dominant energy 
conditions (0 < P < p), Q./2 < jyd?xVhp/T < n. Hence, by Eq. 

1 < ^Hl < 1 (2Q) 

2 - M + 2 (Pyy/^LP/n " ' ^ ' 

For a boundary that corresponds to a stellar surface, P = 0, and Eq. ([26] simplifies, ^ <^ < 
The results above suggest that the Komar mass is a strong candidate for the internal energy in a 
self-gravitating system. However, at the moment this assertion is only tentative: we have found no 
natural splitting of the one-form G into internal energy and work. In any case, the internal energy is not 
an independent variable in the thermodynamic state space of self-gravitating systems — unlike ordinary 
thermodynamics — but a function of the more fundamental geometric variables and Hij . 
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